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ABSTRACT 


The  effect  of  surface  roughness  on  turbulent  heating  rates  is  studied, 
using  a second-order  closure  model  which  specifically  describes  the  effect 
of  roughness  elements  on  the  boundary  Layer.  Comparisons  with  detailed 
measurements  from  low  speed  flat  plate  tests  are  presented  to  verify  the 
model.  A significant  number  of  successful  comparisons  have  also  been  made 
with  wind  tunnel  heat  transfer  data,  mostly  for  hemispherical  noses.  The 
rough  wall  computations  are  analyzed  to  determine  the  nature  of  the  rough- 
ness influence,  and  to  derive  improved  heat  transfer  correlations  for 
engineering  applications. 
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I.  INTRODUCTION 


It  is  well  known  that  surface  roughness  can  cause  significant  increases 
in  turbulent  skin  friction  and  heat  transfer.  An  understanding  of  roughness 
effects  is  essential  to  accurate  design  predictions  for  a wide  variety  of  ap- 
plications, including  ship3,  aircraft,  turbine  blades,  missiles,  and  reentry 
vehicles.  The  present  study  was  motivated  primarily  by  the  last  application. 
At  low  flight  altitudes,  the  thickness  of  the  boundary  iayer  on  the  blunted  nose 
region  of  a hypersonic  reentry  vehicle  can  easily  be  less  than  the  inherent 
surface  roughness  of  practical  heatshieid  materials,  and  roughness  effects 
dominate  the  heat  transfer. 

Most  available  models  for  the  influence  of  surface  roughness  on 

boundary  layer  behavior  are  essentially  extensions  of  Nikuradse' s study  of 

pipes  roughened  with  sand,*  and  the  application  of  Nikuradse' s results  to  flat 

2 3 

plates  by  Prandtl  and  Schlichting.  Dvorak  used  integral  methods  in  which 
the  skin  friction  coefficient  is  specified  as  a function  of  boundary  layer  thick- 
ness and  roughness  height.  Using  this  specification,  moment  equations  are 

4 

solved  for  the  displacement  and  momentum  thicknesses.  Chen  extended 
this  approach  to  predict  heat  transfer,  by  using  the  Stanton  number  correla- 
tion derived  from  subsonic  data  by  Owen  and  Thomson?  In  this  approach, 
the  stagnation  enthalpy  profile  was  equated  to  the  velocity  profile.  A similar 
model  has  recently  been  developed  for  reentry  vehicles  by  Dahm  et  al? 

Here  again,  a momentum  integral  approach  is  used,  with  the  skin  friction 
and  heat  transfer  coefficients  based  on  correlations  of  the  low  speed  flat 
plate  data  of  Healzer  et  alJ  and  the  Mach  2.  9 flat  plate  measurements  of 

Q 

Reda  (cf  only).  The  roughness  augmentation  of  heat  transfer  was  found  to 
be  about  60%  of  the  skin  friction  augmentation,  although  it  was  noted  that 
the  range  of  k/9  values  covered  by  the  data  was  somewhat  limited. 

A somewhat  more  involved  approach  was  taker  by  Saffman  and 

Q 

Wilcox,  utilizing  a two-equation  turbulence  model.  however,  the  effect  of 
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roughness  was  treated  rather  empirically,  by  making  the  boundary  condition 
for  pseudo- vorticity  at  the  wall  a function  of  roughness  height.  This  depen- 
dence was  derived  so  as  to  fit  the  observed  variation  of  the  "law  of  the  wall" 
velocity  deficit  with  roughness.  Some  encouraging  profiles  were  also  com- 
puted for  the  mean  and  fluctuating  velocities.  However,  heat  transfer  was 
again  determined  by  invoking  a Reynolds  analogy  with  the  skin  friction. 

The  present  study  employs  a second-order  turbulence  model  with 
several  attractive  features.  Most  importantly,  it  contains  a fundamental 
description  for  the  effect  of  roughness  elements  on  the  boundary  layer. 

Also,  the  dependent  variables  include  fluctuating  temperatures  as  well  as 
velocity  fluctuations,  so  that  no  assumptions  have  to  be  introduced  relating 
the  heat  transfer  and  skin  friction.  Finally,  solutions  are  obtained  by 
finite-difference  techniques,  integrating  from  the  wall  (y  - 0),  with  no  as- 
sumptions regarding  profiles  or  behavior  in  the  limit  y ~*  0.  Computed  re- 
sults will  be  compared  with  a variety  of  low  speed  and  supersonic  data,  and 
the  nature  of  the  roughness  effect  on  friction  and  heat  transfer  will  be  des- 
cribed. Suggestions  will  be  advanced  for  improved  engineering  specifications. 
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II.  ROUGH  WALL  TURBULENCE  MODEL 


The  turbulence  model  used  here  is  one  in  which  closure  approxima- 
tions are  applied  at  the  second  - or  der . With  the  exception  of  the  treatment 
of  roughness  (to  be  discussed  in  detail  below),  the  formulation  is  somewhat 
standard  at  this  time,  and  has  been  successfully  applied  to  a variety  of 
smooth  wall  boundary  tayer  and  free  shear  flows.  The  model  accounts  for 
both  mean  and  fluctuating  velocities  and  temperatures.  The  dependent  ve- 
locity variables  are  the  mean  velocity  vector  U.,  the  Reynolds  stress  tensor 
u.'u.',  and  the  isotropic  dissipation  rate  5.  Under  the  boundary  layer  ap- 

i | _ _ — 

proximation,  this  set  of  variables  reduces  to  U,  V,  u'  , v'  , w1  , u' v' , and 

2 2 2 

J.  In  practice,  it  is  convenient  to  replace  u'  , v'  , w'  by  the  kinetic  energy 


’222 
l“  (u1  tv'  + w' 


)/2  and  two  measures  of  the  degree  of  anisotropy 


2/3  q , v'  - 2/3  q . 

For  steady  flow,  the  governing  equations  include  continuity: 


the  mean  momentum  equation: 
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and  Re,  is  the  turbulent  Reynolds  number  qA/v,  with  A being  related  to  the 
A 

dissipation  rate  by 
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The  various  terms  on  the  right  side  of  Eqs.  3-7  describe,  respectively, 
the  processes  of  production,  dissipation,  turbulent  diffusion,  molecular  diffu- 
sion, and  streamwise  acceleration.  The  required  closure  approximations 
have  been  developed  from  basic  laboratory  turbulence  experiments,  wherever 
possible  (see  refs.  10-13  for  uetails).  We  shall  not  describe  the  results  here, 
except  to  note  that  they  are  similar  to  those  of  Launder  et  al/^  ^ 

For  flows  with  compressibility  or  heat  transfer,  the  analogous  variables 
involving  the  temperature  (or,  more  precisely,  the  static  enthalpy  h)  are  re- 
quired. These  consist  of  the  static  enthalpy  h,  the  mean  square  fluctuating 
enthalpy  h'  2,  and  the  transverse  and  axial  components  of  the  Reynolds  heat 


flux  v'h'  and  u'h'.  Required  closure  approximations  have  been  made  in  a man- 
ner analogous  to  those  leading  to  the  velocity  equations  above.  The  resulting 
enthalpy  equations  are: 
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It  should  be  noted  that  terms  involving  fluctuating  densities  ( p' ) have 
been  dropped  in  deriving  Eqs.  (1)  - (12).  i’his  is  generally  permissible  if 
the  edge  Mach  number  is  below  4 or  5,  as  is  usually  the  case  for  nosetip 
regions.  However,  the  dominant  effects  of  density  fluctuations  can  be  in- 


cluded by  defining  a generalized  Reynolds  stress  R„  piuu./p 


1 1 


p'  u|  u'.  / p.  Once  this  is  done,  the  primary  effect  of  density  fluctuations  is 
contained  in  a relatively  unimportant  diffusional  term  involving  p'v1,  which 
can  be  related  to  v'  T' . The  resulting  formulation  has  yielded  good  compar- 

17  . 

isons  with  the  measurements  of  Horstman  et  al.  in  a boundary  layer  at 


M 


7. 


The  terms  R , R , R.  contain  the  effect  of  roughness  on  the  bound- 
u q $ 

ary  layer.  Only  distributed  roughness  (comparable  to  sand-grain  roughness) 
is  considered  here,  and  we  make  the  fundamental  assumption  that  the  flow 
around  individual  elements  is  attached  to  the  elements.  For  two-dimensional 
roughness,  or  closely  packed  distributed  roughness,  the  flow  might  be  treated 
more  appropriately  as  cavity  flows  between  the  elements.  (It  must  be 
pointed  out  that  the  present  approximation,  in  which  the  flow  approaching  an 
element  is  assumed  to  be  attached  and  parallel  to  the  wall  even  after  having 
flowed  past  many  upstream  elements,  can  be  improved  upon  for  many  rea- 
listic surfaces.  ) 


i 

i 
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Roughness  elements  provide  a distributed  sink  (due  to  drag)  for  mean 
momentum,  and  distributed  sources  for  kinetic  energy  and  dissipation.  We 
idealize  the  rough  surface  as  being  made  up  of  identical  elements  (although 
the  extension  to  a size  distribution  is  feasible),  with  simple  shapes  such  as 
cones  or  hemispheres.  The  bottom  of  the  elements  corresponds  to  y 0. 

Let  k be  the  element  height,  D(y)  the  element  diameter  at  height  y (for  0 
^ y i k),  and  l be  the  average  center -to-center  element  spacing.  Then, 
viewing  flow  around  an  element  at  height  y as  two  dimensional,  the  form 
drag  between  y - 6y/2  and  y + by/2  is 

~ pU2  CD  D(y )6y 

-2 

To  relate  this  to  drag/unit  volume,  we  note  that  there  are  l elements  per 

2 

unit  surface  area,  so  that  the  appropriate  differential  volume  is  l by  and 
the  sink  term  for  mean  momentum  is 

R ^ pU2  C DU2  (13) 

u 2 D 

We  could  specify  1 for  the  drag  coefficient,  appropriate  to  infinite 

circular  cylinders  at  local  Reynolds  numbers  above  the  Stokes  flow  regime. 
However,  lower  values  such  as  0.  5 are  more  appropriate  for  finite  elements 
such  as  cones  and  hemispheres,  and  use  of  such  a value  provides  a first- 
order  correction  for  the  three-dimensional  nature  of  the  flow  about  the 
element  s. 

Specification  of  the  source  terms  R and  R is  more  speculative,  but 

q $ 

these  terms  are  generally  smaller  than  the  natural  production  terms  and 
hence  less  critical  to  the  model.  Fluctuations  are  introduced  as  disturbances 
in  the  wakes  of  the  elements.  By  analogy  with  familiar  wake  flows,  fluctuat- 
ing velocity  components  should  be  on  the  order  of  two-tenths  of  the  local  flow 
velocity  (see,  for  example,  Ref.  18).  Thus,  the  kinetic  energy  created  per 
unit  volume  is 


'I 


R 0.  04  pU3  D/22 

q 


U) 


A low  Reynolds  number  cutoff  to  this  term  was  also  considered,  anticipating 
that  the  wakes  of  elements  should  be  laminar  below  some  minimum  value  of 
pUD/p.  However,  in  such  situations,  the  local  flow  velocity  II  is  sufficiently 
small  that  R is  already  negligible. 

1'he  role  of  the  source  term  for  dissipation  R,  is  to  control  the  wave- 

J 

length  of  the  fluctuations  created  by  the  roughness  elements.  Again,  by  the 
wake  analogy,  this  wavelength  should  be  comparable  to  the  diameter  of  the 
element.  Then,  following  the  reasoning  leading  to  Eq.  14, 


Rx  0.  04  pU3v/(D22) 
$ 


(15) 


In  evaluating  the  calculations  presented  below,  this  term  was  found  to  be 
negligible.  However,  the  term  is  given  here  for  the  sake  of  completeness. 

1’he  distributed  source  or  sink  terms  given  by  Eqs.  (13),  (14),  and 

(15)  are  the  only  ones  that  need  be  considered.  If  oscillations  in  the  wake 

of  an  element  are  approximately  isotropic,  there  should  be  no  significant 

2 2 J 

creation  of  the  anisotropic  components  u'v’,  Sj  ^ u1  - 2/3q  , S v1 
- 2/  3 q2.  Except  in  the  Stokes  flow  regime,  heat  transfer  to  an  element 
should  be  small.  Therefore,  there  should  be  no  distributed  source  or  sink 
terms  in  the  equations  for  the  thermal  variables. 

Boundary  conditions  to  Eqs.  (1)  - (12)  are  obvious:  fluctuating 
quantities  are  zero  at  a solid  wall  or  at  the  outer  edge  (if  there  is  no  free- 
stream  turbulence).  For  numerical  solutions,  the  equations  are  first 
transformed  to  the  standard  str eamfunction  coordinate,  guaranteeing  con- 
tinuity and  eliminating  the  normal  velocity  V.  The  transverse  coordinate 
is  also  normalized  by  the  edge  value  of  the  stream  function,  so  that  addi- 
tional mesh  points  need  not  be  carried  in  the  free  stream  to  allow  for 
boundary  layer  growth.  For  proper  resolution  of  the  region  near  the  wall, 
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a linear  mesh  in  the  logarithm  of  the  str eamfunction  is  used.  The  finite- 
difference  equations  are  solved  with  a block  tridiagonal  Newton-Raphson 
technique. 


III.  INCOMPRESSIBLE  TURBULENT  BOUNDARY  LAYERS 


The  basic  case  of  a low  speed  boundary  layer  on  a flat  plate  is  im- 
portant to  examine,  if  only  because  detailed  measurements  are  available. 

7 19  20 

In  particular,  the  experiments  of  Moffat  and  co-workers  ’ ’ on  the  flow 

over  a plate  covered  with  0.  05  inch  diameter  spherical  balls  provide  data  on 
skin  friction,  heat  transfer,  and  profiles  of  the  mean  and  fluctuating 
quantities. 

Figure  1 compares  the  results  of  our  model  with  the  mean  velocity 
data,  plotted  in  law-of-the-wall  coordinates.  Obviously,  there  is  a signifi- 
cant shift  in  the  mean  velocity,  although  the  smooth  wall  and  rough  wall 
profiles  are  more  alike  when  plotted  against  y/6.  The  three  components  of 
the  kinetic  energy  are  «hown  in  Fig.  2,  and  Fig.  3 gives  the  Reynolds 
stress  profile.  Perhaps  the  only  discrepancy  occurs  for  the  rms  velocity 
component  in  the  flow  direction  (u' ) near  the  wall.  However,  this  effect 
seems  to  occur  with  all  second-order  models,  even  on  smooth  walls.  It 
has  yet  to  be  explained  properly,  but  seems  not  to  affect  any  other  boundary 
layer  properties  of  interest.  The  roughness  effect  is  apparently  substantial 
in  Figs.  2 and  3;  much  if  this  could  be  scaled  out  by  normalizing  the  fluctuat- 
ing velocities  by  the  friction  velocity  rather  than  the  edge  velocity. 

Finally,  the  skin  friction  coefficient  and  Stanton  number, 

q/p  u (H  - H ),  are  shown  as  a function  of  downstream  distance  and  flow 
e e e w 

velocity  in  Figs.  4 and  5.  Increasing  the  flow  velocity  leads  to  a thinner 
boundary  layer  and  a larger  effect  of  roughness.  Although  the  smooth  wall 
case  was  not  investigated  by  Healzer  et  al J , the  present  calculations  are  in 
good  agreement  with  classical  values.  Again,  the  model  agrees  veil  with 
the  data,  and  it  is  noted  that  the  effect  of  roughness  on  skin  fricticn  is 
greater  than  that  on  heat  transfer.  For  the  case  with  the  lowest  velocity 
(32  ft/sec),  transition  was  apparently  not  far  upstream  of  the  measured 
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Fig.  1 Flat-Plate  Velocity  Profiles  for  Doth  Rough 
and  Smooth  Walls. 
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DIMENSIONLESS  DISTANCE  FROM  THE  WALL 


VELOCITY  FLUCTUATION  PROFILES  » 10 
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Fig.  2 Velocity  Fluctuation  Components  for  Both 
Rough  and  Smooth  Walls. 
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Fig.  3 Profiles  of  the  Reynolds  Stress  for  Both 
Rough  and  Smooth  Walls. 
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Fig.  4 Skin  Friction  Coefficient  on  a Roughened 
Flat-Plate  as  Function  of  Downstream 
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Fig.  5 


Normalized  Heat  Transfer  (Stanton  Number) 
Distribution  as  Function  of  Downstream  Distance 
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IV.  SUPERSONIC  BLUNT  NOSE  HEAT  TRANSFER 


Of  more  practical  interest  is  the  heat  transfer  to  the  blunted  nose 
region  of  a high  speed  flight  vehicle.  At  sufficiently  high  Reynolds  numbers, 
corresponding  to  low  flight  altitudes,  the  transition  location  can  be  quite 
close  to  the  stagnation  point.  Also,  the  boundary  layer  momentum  thickness 
can  be  as  small  as  one  mil,  so  that  only  a highly  polished  surface  would  be- 
have in  a smooth  manner. 

The  data  considered  in  this  section  were  obtained  by  the  Acurex  Corp. 

2 1 

under  the  PANT  program,  in  NSW C Tunnel  No.  8 at  a freestream  Mach 
number  of  5.  The  thin-wnll  calorimeter  models  generally  had  a nose  radius 
of  2.  5 in.  Various  degrees  of  surface  roughness  were  generated  either  by 
grit  blasting  or  by  brazing  particles  onto  the  surface;  the  resulting  roughness 
is  thought  to  be  similar  in  character  to  sand  grain  roughness,  and  the  quoted 
roughness  heights  are  "peak  to  valley"  values  (which  might  be  slightly  larger 
than  the  equivalent  sand  grain  roughness).  The  ratio  of  wall  temperature  to 
freestream  total  temperature  was  in  the  range  0.4-  0.  7. 

Figures  6 through  9 show  comparisons  with  the  PANT  heat  transfer 
coefficient,  q/(T  - T ),  data  ver  sus  di stance  from  the  stagnation  point, 

6 W 

for  decreasing  freestream  unit  Reynolds  number.  The  computations  were 
started  near  the  stagnation  point  (S/R^  < 0.  01)  with  a laminar  profile.  The 
effect  of  roughness  generates  a turbulent  boundary  layer  relatively  quickly, 
although  this  procedure  should  not  be  considered  an  accurate  prediction  of 
transition.  In  some  cases,  such  as  the  smooth  case  of  Fig.  6,  the  local 
Reynolds  number  is  insufficient  for  establishment  of  a turbulent  boundary  layer 
for  some  distance  from  the  stagnation  point  (S/R^  ~ 0.  25  in  that  case). 

It  appears  from  Figs.  6 through  9 that  the  theory  generally  provides 
for  good  agreement  with  the  measurements.  A few  points  should  be  made, 
however:  i)  the  smooth  wall  case  was  almost  certainly  not  fully  turbulent, 
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Fig.  7 


Comparison  With  PA 
Re  = 10  x 106 /ft. 
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HEAT  TRANSFER  COEFFICIENT  (BTU/ft  -sec 


r-ven  at  the  highest  Reynolds  number;  ii)  at  larger  angLes  (S/R^  > 1),  the 

computations  are  inaccurate  due  to  breakdown  of  the  Newtonian  pressure 

distribution  employed;  iii)  significant  apparent  increases  over  the  theoretical 

laminar  heating  rate  are  observed  at  the  stagnation  point  (S/R^  = 0)*  The 

approach  taken  here  yields  the  laminar  result  at  the  stagnation  point  where 

there  is  no  production  of  turbulence  since  U =0  and  5U/5Y  = 0.  We  in- 

e 

vcstigated  whether  a low  Reynolds  number  (Stokes)  drag  coefficient  in  R 
(Eq.  (13))  plus  the  corresponding  heat  transfer  term,  could  increase  the 
heat  transfer  near  (but  not  at)  the  stagnation  point.  The  result  was  negative. 
In  our  opinion,  the  effect  of  roughness  on  the  stagnation  point  heat  transfer 
is  not  presently  understood.  Random  motion  of  the  stagnation  point,  due 
perhaps  to  tunnel  disturbances,  would  increase  the  heating  rate  at  the  mean 
stagnation  point  location.  However,  for  the  effect  to  be  appreciable,  the 
amplitude  of  such  motion  would  have  to  be  ~ 10  - 20°,  which  seems  unlikely. 
The  stagnation  point  could  be  unstable  to  3-D  disturbances  of  the  Goertler 
type,  but  the  role  of  roughness  in  triggering  such  instabilities  has  not  been 
demonstrated.  We  suspect  that  the  stagnation  point  data  of  Figs.  6 through 
9 are  contaminated  by  heat  conduction  within  the  model.  A lateral  conduction 
correction  was  applied  in  the  PANT  data  reduction,  although  the  thermocouple 
spacing  is  much  too  coarse  to  resolve  the  issue. 
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V.  HEAT  TRANSFER  SCALING  LAWS 


One  particularly  useful  aspect  of  this  rough  wall  turbulence  model  is 
that  the  results  can  bo  examined  to  determine  the  nature  of  the  roughness 
influence  on  turbulent  boundary  layers.  One  rather  conspicuous  conclusion 
is  that  the  Reynolds  analogy  between  friction  and  heat  transfer  is  not  pre- 
served with  significant  roughness.  This  result  is  well  known,  and  derives 
from  the  absence  of  a heat  transfer  analogy  to  form  drag  on  the  elements. 
The  computations  show  that  the  velocity  fluctuations  increase  in  pr oportion 

to  II  ( \ T / p ) , but  the  temperature  fluctuations  are  hardly  changed  by 

T \V  ® 

roughness.  Since  T ~ u'v1  and  q~  v'  T'  the  heat  Lransfer  augmentation  is 

w 

the  square  root  of  the  skin  friction  augmentation: 


o 


(lh) 


where  subscript  o denotes  smooth  wall. 

To  indicate  the  appropriate  scaling  law  for  the  roughness  influence, 
the  wall  shear  is  formally  given  by 


T 

W 


w 


r 2 i 

pU  D(y ) — dy 

Jo  iZ 


(17) 


As  indicated  by  Fig.  1,  the  velocity  is  significantly  altered  by  roughness. 
However,  for  the  present  purposes  we  will  use  the  smooth  wall  law-of-the 
wall  to  evaluate  the  integral  in  Eq.  (17): 

U U (2.  5 in  y+  t 5.  1)  (18) 

To 

In  addition  to  being  in  error  by  up  to  a factor  of  two  due  to  roughness,  the 
roughness  height  can  easily  extend  into  the  wake  region,  beyond  the  validity 
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ol'this  relation.  Neglecting  contributions  from  the  laminar  sublayer,  the 

lower  limit  of  the  integral  should  be  replaced  by  y 10.  For  simplicity, 

we  set  P(y)  k and  p p in  Fq.  (17).  Inserting  Eq.  (18)  into  Eq.  (17) 

\v 

and  rearranging. 


T t:  p 

— 1 I — f.(k‘) 

T t.  pl 

"o  fo 


f.(k')  ^*7-7  fk  • or,!,1 

**  X •*  k '10 


(19) 


4^  |in. 


t 2*1 

k I 6.  25  Xn  k - t,3 


(20) 


i 

where  k p F k/p  . 

w ~ w 

Figure  10  shows  the  roughness  augmentation  of  skin  friction  in  the 
form.it  of  Eq.  (19),  according  to  our  computations.  Agreement  with  the 
incompressible  flat  plate  data  was  sufficient  (cf.  Fig.  4)  to  ensure  that  the 
data  would  plot  in  a similar  manner,  but  the  skin  friction  was  not  measured 
in  the  sphere  tests.  The  plotted  points  correspond  to  fully  turbulent  condi- 
tions; values  for  S/R  < 0.  2 or  Re.  < 150  were  omitted.  Fhe  solid  line  on 

IN  o 

the  figure  is  Eq.  (20),  with  0.  5 and  X 5k.  Given  the  assumptions  re- 

quired to  reach  Eq.  (20),  the  agreement  is  surprisingly  good.  However,  a 

wider  range  of  data  should  be  examined  before  Eq.  (20)  would  be  considered 

2 2 

reliable.  In  particular,  the  roughness  density  effect  given  by  k /X  is  prob- 
ably oversimplified.  Changes  in  roughness  density  would  be  expected  to 
yield  variations  in  the  actual  velocity  profile  and  could  lead  to  a more  com- 
plex dependence  on  roughness  density. 

Figure  11  shows  the  corresponding  plot  for  the  heat  transfer  aug- 
mentation, again  based  entirely  on  our  computed  values.  The  line 
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Fig.  10  Roughness  Augmentation  of  Skin  Friction 
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Fig.  11  Roughness  Augmentation  of  Heat  Transfer 
Open  symbols:  k/t'<  1: 

Filled  symbols:  k/0>  l. 
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correspond.-  to  t f 1^  ; this  manner  of  applying  the  square  root  is  only  an 

approximation  to  Kq.  (lt>),  and  may  be  partially  responsible  for  the  greater 
scatter  in  '.his  figure  compared  to  the  previous  one.  Another  reason  for  the 
scatter  is  that  Kq.  (lt>)  is  not  strictly  correct,  in  that  the  temperature 
fluctuations  can  be  affected  somewhat  by  roughness. 

By  way  of  comparison,  the  skin  friction  expressions  used  by  Dvorak 
and  Chen  are  expressed  in  terms  of  the  displacement  thickness  rather  than 
wall  coordinates: 

~ A <n6  /k  i f(density)  . 


i'lie  dependence  on  roughness  density  has  yet  to  be  evaluated  with  our  model; 
much  of  the  data  from  which  their  expression  was  derived  pertain  to  two- 
dimensional  roughness,  which  was  not  considered  here. 

The  correlations  of  Dahm  et  al^  are  closer  to  the  present  result, 
although  still  with  important  differences.  Their  expressions  are: 


0 


1 ! 0.  5 f(k/e  )g(  x) 


(21) 


1 t 0.  3 f(k / 9 )g( x) 

bto 

k / 

f ( k / 0 ) l t 0.  09  k/0  + 0.  53  (1  - e ) 

t 

g(x)  x t 1.5  (1  - e X)  , x log  jq  Ys  5 


(22) 


The 

ing  since  it 


function  g(x)  is  quite  similar  to  our  f^(k  ),  which  is  not  surpris- 
represents  a correlation  of  essentially  the  same  data.  The 
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dependence  on  k/0  is  weak:  f(k/0)  varies  by  only  about  70%  between 
k/t?  1 and  10.  We  did  not  consider  a dependence  on  k/6,  although  it  might 
arise  at  large  roughness  heights  that  extend  beyond  the  logarithmic  velocity 
region.  However,  as  indicated  by  the  open  and  filled  symbols  in  Fig.  10, 
no  strong  effect  is  evident.  The  difference  between  the  factors  of  0.  5 and 
0.  t in  Fq.  ( 2 1 ) and  (22)  accomplishes  much  the  same  effect  as  the  square 
root  in  our  result,  Eq.  (16).  One  difference  that  is  somewhat  significant  is 
the  lack  of  a density  dependence  in  the  Dahm  correlation;  this  dependence 


shows  up  rather  clearly  between  the  subsonic  flat  plate  cases  (p  /p 

w e 

/v 

and  the  supersonic  sphere  cases  (p  /p  0.  5). 

w e 
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VI.  FINAL  REMARKS 


The  results  presented  here  are  quite  encouraging,  in  that  a rather 
basic  model  yields  results  that  are  in  agreement  with  many  observed  trends 
regarding  the  influence  of  surface  roughness.  l’he  assumptions  inherent  to 
this  m >del  are  limited  to  the  basic  nature  of  the  flow  around  the  roughness 
elements,  and  no  approximations  have  been  made  regarding  profiles  of  the 
boundary  layer  quantities,  turbulence  levels,  or  relations  between  the 
momentum  and  energy  fluxes. 

Flu*  present  results  can  undoubtedly  be  improved  and  extended.  It 
would  be  useful  to  have  better  data  to  study  on  the  effects  of  roughness  den- 
sity and  element  shape,  for  distributed  roughness.  Also,  the  scaling  laws 
derived  in  the  previous  section  could  be  improved  if  we  could  better  correlate 
the  mctn  velocity  profiles  in  the  presence  of  roughness.  Finally  there  are 
other  situations  that  could  be  examined  with  the  approach  presented  here, 
including  the  combined  effect  of  roughness  and  mass  addition,  and  the  effect 
of  roughness  at  strongly  supersonic  or  hypersonic  edge  Mach  numbers. 
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